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Boltzmann-Gibbs distribution arises as the statistical equilibrium probability distribu- 
tion of money among the agents of a closed economic system where random and undi- 
rected exchanges are allowed. When considering a model with uniform savings in the 
exchanges, the final distribution is close to the gamma family. In this work, we imple- 
ment these exchange rules on networks and we find that these stationary probability 
distributions are robust and they are not affected by the topology of the underlying 
network. We introduce a new family of interactions: random but directed ones. In this 
case, it is found the topology to be determinant and the mean money per economic agent 
is related to the degree of the node representing the agent in the network. The relation 
between the mean money per economic agent and its degree is shown to be linear. 
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Introduction 

Agent-based modeling can be used to study systems exhibiting emergent prop- 
erties which cannot be explained by aggregating the properties of the system's 
components. 1 Statistical mechanics and economics share the property to analyze 
big ensembles where the collective behaviour is found out as a result of interac- 
tions at the microscopic level and where agent-based simulations can be applied. 
Many systems are studied in terms of the nature that defines their inner compo- 
nents while others are considered from the point of view of the interactions among 
the agents that can be pictured through a complex network. Plenty of informa- 
tion is encoded in connectivity patterns. Hierarchical structures appear in a natural 
way when we study societies and, according to many authors^ one of the mile- 
stones is to understand why and how from individuals with initial identical status, 
inequalities emerge. This is related to the question of hierarchy formation as a self- 
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organization phenomenon due to social dynamics: - Elitarian distributions can arise 
starting from a society where people initially own an equal share of economic re- 
sources, e.g.: the exponential distribution for the low and medium income classes 
in western societies! 4 ! 5 ! 6 ! 

We consider Dragulescu-Yakovenko gas-like models in economic exchanges^ so, 
let our system be composed of ./V economic agents, being N 3> 1 and constant. Each 
agent i owns an amount of money m, so the state of the system at a given time is 
defined by the values that every variable rrii takes at that moment, {rrii}^. Money 
distribution among the agents should never be confused with the notion of wealth 
distribution. Money is only one part inside the whole concept of wealth. Transfer of 
money represents payment for goods and services in a market economy. We study 
simplified models which keep track of that money flux but do not keep track of what 
goods or services are delivered. At each interacting step, agents trade by pairs and 
local conservation of money is sustained, 

(mi, rrij) i — > (to-, m'j) : m\ = rrn + Am, (1) 

fflfj = rrij — Am. 

Transactions result in some part of the money involved in the interaction changing 
its owner. For simplicity, we do not consider models where debts are allowed. 

It is deeply established in the common knowledge that highly-ranked individuals 
in societies have easier access to resources and better chances to compete. This is 
a motivation to look for internal correlations between money and surrounding en- 
vironment. We wonder if the exchange rules that define simple gas- like models for 
random markets, when implemented on networks, are capable of depicting correla- 
tions between purchasing power of an agent inside a social network and the influence 
of the agent on the rest of the system. We associate the purchasing power concept 
to the mean money per economic agent computed as a function of the connectivity 
degree of each agent in a network. At this level, influence of an agent is only related 
to the degree of the node representing the agent. We implement the exchange rules 
on two type of networks: uniform random spatial graph and Barabasi- Albert model, 
and then examine the relationship between the former econo-social agent indicators 
for the different underlying architectures. 

In section [TJ we review two well-known random undirected exchange rules: gen- 
eral and uniform savings models. In section [2j we introduce a new family of inter- 
actions: random but directed ones. The main property of this simple exchange rule 
is that it is a real inspired model where social inequalities in money distribution 
emerge in a natural way. In section [3j we show the relation between mean money 
per economic agent and the connectivity degree of the agent. For the models with 
undirected exchange rules, we observe no correlation between money and the degree 
of the nodes. Linear dependence is found for the new random exchange model we 
propose. Section [4] is devoted to gather the most relevant conclusions. 
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1. Undirected random markets 
1.1. Undirected random market 

For some random economic systems where money is a conserved quantity, the 
asymptotic distribution of money among the agents is given by the Boltzmann- 
Gibbs distribution (BG), 

p cq ( mi = m ) = -L e -™/(™), (2 ) 
(to) 

where the role of the effective temperature is played by the average amount of money 
per agent, 

1 N 

i=l 

This feature was first shown by Dragulescu and Yakovenko in 2000 by means of nu- 
merical simulationsP Subsequently, analytical justification was given by Lopez- Ruiz 
et al. in 2008 and 2012. BG can be geometrically deducecP under the assumption of 
equiprobability of the possible economic microstates. We also know that an asymp- 
totic evolution towards BG is obtained regardless of the initial distribution for those 
systems with total money fixed and when considering random symmetric interac- 
tions between pairs of components.^ This comes from BG being the stable fixed 
point of the distributions' space, Lf[0, oo) = {p(x) : J °° p(x) dx < oo}, under the 
iterated action p(x) — > p'(x) of the integral operator T given by 

p'( z )=[Tp](z)= [[ P^ldxdy, (4) 

JJs(z) x + y 

where S(z) = {(x,y), x,y > 0, x + y > z} is the integration domain. 

Let us now consider the gas-like model originally proposed by Dragulescu and 
Yakovenko so, at each computational step, we randomly choose a pair of agents and 
then, one -the labeled as i- is chosen to be the winner in the interaction process 
and the other one -labeled as j- becomes the loser and, according to the previously 
stated rule ([l]), an amount of money Am is transferred from the loser to the winner. 
Assuming Am > 0, it is obvious that if the loser does not have enough money to 
pay, which is nothing but the local condition rrij < Am, the transaction is forbidden 
and we should proceed with a different pair of agents. 

Instead of considering the restriction in the interaction, we state the exchange 
rule considering Am — e(TOj + mj) — m, that gives rise to the completely random 
case given by 

(m, mj) ' — > (m'i, m'j) : (5) 
m- = £(m, + mj), 
m'j = (1 — e)(m, L + m,j), 
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Fig. 1. Undirected random market. Numerical simulation considering N = 5000 agents and 
averaging over 250 samples. Every agent starts with mj = (m) = 0.5 and the exchange rule is given 
by equation |HJ. Stationary probability distribution of money P(m), computed after N 2 interac- 
tions, compared to the blue solid curve describing the Boltzmann-Gibbs law P(m) <x e~ m l^ m \ 
Log-linear plot and entropy time-evolution are the insets. 

where s € [0, 1] is a uniform random number which is refreshed at every computa- 
tional step. Observe that both agents can be winner or loser in a symmetric way, 
depending on the random number e at each step. This approach also ensures that no 
agent will evolve to own a negative amount of money or, in other words, debts are 
not allowed. The condition m, > for every agent i in the system is accomplished in 
a natural way. This exchange rule ^ is a very rough macroeconomic model where 
individuals or corporations raise their money for a venture and then, the market 
effect or their mutual interaction determines the final distribution. See figure [I] 

1.2. Undirected random market with uniform savings 

The concept of savings arises in an obvious way from observing human behaviour 
when people are inmerse in a market economy.^ This feature is introduced through 
a parameter, A € [0,1], which is called a propensity factor .121 This means that 
each agent saves a fraction A of its money when an interaction occurs and trades 
randomly with the other part: 

(rrii, rrij) \ — > (m*, m'j) : (6) 
m! i = \rrii + e(l — A)(m.j + rrij), 
m'j = Xnij + (1 — e)(l — A)(mj + n%j). 

We consider the model with uniform savings which means that A is fixed to be 
constant among the agents and with no dependence on the time. The statistically 
stationary distribution P(m) decays rapidly on both sides of the most probable 
value for the money per agent which, in this case, is shifted from the poorest part of 
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Fig. 2. Undirected random market with uniform savings. Numerical simulation considering 
N = 5000 agents and averaging over 250 samples. Every agent starts with m; = (m) = 0.5 and 
the exchange rule is given by equation JHJ. Stationary probability distribution of money P(m), 
computed after N 2 interactions for different values of A. Entropy time-evolution insets. 

the system to (to) when A — > 1. See figure [2j This behaviour was already described 
as a self-organising feature of the market induced by self-interest of saving by each 
agent without any global perspective in an analogous way to the self-organisation 
in markets with restricted comoditiesP^ 

First attempt to give a quantitative description for the steady distribution to- 
wards this model evolves is due to Patriarca et al. in 2004. They stated that nu- 
merical simulations of ([6j could be fitted to a standard Gamma distribution!^ 
Subsequently (2007), Chatterjee and Chakrabarti offered a brief study of the con- 
sequences that this modelization implies and stated that as A increases, effectively 
the agents retain more of its money in any trading, which can be taken as implying 
that with increasing A, temperature of the scattering process changes.^ According 
to their study, fourth and higher order moments of the distributions are in discrep- 
ancy with those of the Gamma family so, the actual form of the distribution for this 
model still remains to be found out. Calbet et al. (2011) gave an iterative recipe to 
derive an analytical expression solving an integral equation.^ A similar expression 
was derived in a different way by Lallouache et al. (2010) ff3 



2. New scenario: directed random market 

We have shown how the propensity factor A is introduced |B| as a variation for 
the general random undirected exchange rule ([5]) showing how, from individual 
responsible decissions -such as saving a fraction of your money when entering an 
exchange market-, self-organized distributions where the mean and the mode are 
close arise so, richness would be quite balanced distributed among the group. A 
completely different scheme^ that modifies the general rule ([l]) proposes a random 
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sharing of an amount 2m j (instead of m, + m,j) only when rrii > m,j, trading at 
the level of the lowest economic class in the trade. This model leads to an extreme 
situation in which all the money in the market drifts to one agent and the rest 
become truely pauper. From this idea, we give a new and more general exchange rule 
reflecting this directed or biased orientation for the interaction and including this 
particular result. We propose an integral operator which is the analytical approach 
to this new rule in the mean-field or gas- like case and, in section [3j we implement 
this rule on networks to study how it is affected when we mix directed interactions 
with undirected networks. 

The directed exchange can be understood as a first approach to microeconomic 
activities where money is transferred only in one direction, similar to payments for 
goods. We consider the most general family of interations, 

(m», Tfij ) i — > (m'i, m/j) : (7) 
= srrii, 
m'j = m,j + (1 — e)mj. 

where e G [0, 1] is a random number chosen with uniform probability. At each time, 
the system is described by the probability distribution function of money when we 
choose one of the agents randomly, P(m). In the continuous limit, we can picture 
the system to be the combination of two identical copies, Pi(u) and P2(v), of the 
original system itself, 

P(m) = ~P x (m) + \P2{m). (8) 

For each interaction we can consider that the two different agents i and j, with 
money values being rrij and rrij, are two realizations of picking up randomly one 
agent from each copy, P\{u = m.j) and /^(^ — m j) respectively. In every interaction 
step (u, v) — > (vf , v 1 ), there is a transaction of the agents conforming the distribution 
P(m) to a new configuration given by 

P\rn) = \p[{m) + \p^m), (9) 

where the agents u' will conform P[(m) and v' will conform P^m). Let us now 
consider the probability of a randomly chosen agent among the first copy P[(u') 
owning an amount of money vf — x after the interaction happened. From Q, it 
is clear that u > x, and as the result v! is uniformly distributed in [0, u] so, the 
probability of obtaining a certain value x is given by X/u. The interaction of pairs 
(u, v) in the first configuration of the system gives rise to the evolution of Pi(u) to 
the following probability P[ of obtaining v! = x: 

p{(x)= ff p(u)p(v) dudv= r p{v)dv x r P M du= f pM du . 

JJu>x u JO Ju>x u Ju>x u 

(10) 

For the second copy, P^v'), we should consider that money of the second agent 
after the interaction, v' — x, should be a value between the initial money it has, v, 
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and the maximum possible amount of money it can have after the interaction, which 
will be associated to get all the money from the first agent so, v + u. Again, the 
length of the segment + is u so, the probability of v' is uniformly distributed 
in that segment and so, the probability to have v' = x will be l/u. The expression 
for the probability P' 2 of having v' = x in the second copy of the system results in 

p(u)p(v) 



du dv. 



(11) 



v<x<u+v 



We explicitly compute the norm, 



f P[(x)dx= [ dx f ^-du^ f du f ^-dx= [ p{u)du = l, (12) 

JO Jo Ju>x u Jo Jo u Jo 

y* oo />oo p 

I P^(x)dx = dx 

Jo Jo Jv<x<u+v 



p(u)p(v) 



du dv 



du 



dv 



u+v 



p(u)p(v) 



and the expected value, 



x P[(x) dx = 



x P^x) dx 



x dx 
up(u) du 
x dx 
du 



p(u) 



du = 



dx = 



du 



p{u) du 
p(u) 



p(v) dv = 1, 



dx 



(13) 



du 



v<x<u-\-v 

/•U+V 

dv 

J V 
1 U 2 



P{u)p{v) 
u 

p(u)p(v) 



du dv 
du dv 



2uv 



du 



p{u) p(v) dv 



up(u)p(v)dv+ / du 2vp(u)p{v)dv 



From results (12 1 and (13), together with definition ([9]), the norm and mean 
value of the distributions are conserved when we define the operator for the directed 
random interaction by 



[Tp]{x) = 



p(u) 



du 



p{u)p(v) 



du dv. 



(14) 



v<x<u+v 



Although we cannot give a proof for the infinite iteration of this operator, we see 
that it piles up the distribution at the lower values for m which gives rise to a very 
impoverished population and a very slightest fraction of too opulent agents. Even 
with such a very narrow initial distribution we choose, the effect of this operator is 
very strong in only a couple of iterations. This result is in perfect agreement with 
the first model of biased interaction^ we mentioned in our previous section. Other 
models with asymetric rules that establish a transition between Boltzmann-Gibbs 
and Pareto distributions can be found in the literature.^ See figure [3j 
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Fig. 3. LEFT: First and second application of the integral operator jl4[ l over a Gaussian distri- 
bution. Dotted curve is the initial Gau ssian distribution N((m) = 0.5, cr = 0.05). The distribution 
resulting from the first application of | |l4| , [Tp](:r), is plotted in dashed lines. Second iteration, 
[T[7~p]] (x), is plotted in solid line. Due to the scale proportions, inset figure shows a detailed com- 
parison between initial distribution and the resulting from the first iteration. RIGHT: Directed 
random market implemented as gas-like model. Numerical simulation considering N = 5000 agents 
and averaging over 250 samples. Every agent starts with = (m) = 0.5 and the exchange rule 
is given by equation (J7J. Stationary probability distribution of money P(m), computed after N 2 
interactions and inset showing the whole entropy-evolution. We can appreciate how the system 
starts evolving towards the state of maximum entropy and then it starts a self-organising process 
which is very slow. When considering extremely long runs of the simulation we observe how the 
entropy falls towards states of negligible entropy, as we should expect from a system where lack of 
information in P(m) is just the permutation of the agent owning all the money. 



3. Simulations on networks 

We stated in the introduction of this work the desire of finding a first model which, 
when implemented on networks, is able to show a relation between money and 
influence of an agent on the system. For this purpose, we simulate the two models 
we have already studied in section [T] and the new one we have presented in section 
[5J We choose two representative cases: the random uniform spatial network (SP) 
and the Barabasi- Albert model (BA). SP is an easy way to build networks with 
Poissonian degree distributions^ and BA is an algorithm for generating random 
scale-free networks following the preferential attachment prescriptionPH 

At the beginning, every agent is given the same amount of money so, the initial 
distribution of money among the agents is written 

f initial (m) =5{m- (m)), (15) 

and we obtain the steady state distribution for the exchange rules ^ and ^ 
implemented on these two different topologies. From the histogram related to the 
distribution of money among the agents, we can consider the entropy associated to 



Directed Random Markets: Connectivity determines Money 9 



that distribution in a discrete form, 

S=-^P A (i)logP A (i), (16) 
i=i 

where A just recalls the coarse graining when computig P(m) as a discrete his- 
togram. For simplicity, we measure S/S max , the entropy as a fraction of the maxi- 
mum which is given for the exponential distribution. 

If the money distribution is affected by the topology of the underlying network, it 
should show some kind of dependence on the degree of the agents. For our purposes 
it is enough for the reader to associate the degree of a node to its number of 
neighbouring agents it can interact with. As we will only use for now undirected 
and static networks, the degree of an agent will be a distinctive feature. We define 
the mean money per economic agent as a function of k, (m)(fc), given by 

, n 2 N k 

(™)(fc) = ^rEE m *( T ) Vi : k ^ k - ( 17 ) 

k T=l »=1 

Note that we compute the mean money of the nodes inside each class of connectivity 
at every step of our simulation and then, we consider the time-averaged mean money 
per node according to the different degrees. Nk denotes how many nodes have degree 
k and our simulations run for N 2 steps. There is no need to be worried about the 
transitory regime disturbing the computed results for (to) (k) because it is negligible: 
less than 0.5% of total computed steps in the worst case (see figure [6]). We also 
compute the standard deviation, a m {k), given by 



m = 



N 2 N k 



\ ^EEhM-ww] 2 as) 

\ fc r=l i=l 



Vt 



In figures mm aad [6] we show, for the rules ^ and Q, respectively, the 
stationary probability distribution P(m) with entropy-evolution of S/S max for the 
whole simulation with inset detail of the transitory regime when it is required. When 
we plot the distributions (m)(fc) and cr m (fc), it is also shown the characteristic degree 
distribution of the networks we implement. 

We clearly see in these figures how rules ^ and ^ are transparent to the 
underlying topology, provoking always the decay of the economic system to the BG 
or Gamma- like distributions, respectively, indistinctly of the different connectivities 
of the agents. We can say that for these types of interactions the economic classes 
are blind respect to the social influence of the agents. By contrary, the new directed 
rule ([7]) separates the agents in economic classes correlated with their connectivities, 
in such a way that more connected agents show a bigger propensity to accumulate 
more money, in this case with a linear relationship between money and connectivity. 
This is a characteristic that, in general and independently of the political system 
installed in the power, seems to be more likely to be found in the reality. 
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(a) Spatial network: 
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(b) Barabasi-Albert model: 





Fig. 4. Undirected random market implemented on networks. Numerical simulation considering 
N = 5000 agents and averaging over 250 networks for each type (a) and (b). Every agent starts 
with mi = (m) = 0.5 and the exchange rule is given by equation |5j. LEFT: Stationary prob- 
ability distribution of money P(m), computed after N 2 interactions, compared to the blue solid 
curve describing the Boltzmann-Gibbs law P(m) oc e~ m l ( m '> and entropy-evolution for the whole 
time of the simulation and inserted detail for the initial transitory. RIGHT: Mean money per eco- 
nomic agent computed according to |l7j and its standard deviation, given by \18\ . Typical degree 
distribution is also shown. 
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(b) Barabasi-Albert model: 
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Fig. 5. Undirected random market with uniform savings implemented on networks. Numerical 
simulation considering N = 5000 agents and averaging over 250 networks for each type (a) and (b). 
Every agent starts with m; = (m) = 0.5 and the exchange rule is given by equation JHf. LEFT: 
Stationary probability distribution of money P(m), computed after TV 2 interactions for different 
values of A and entropy-evolution for the whole time of the simulation. Transitory towards equilib- 
rium entropy is analogous to the one shown in figure [4] RIGHT: Mean money per economic agent 
computed according to \17\ and its standard deviation, given by | |18} . Typical degree distribution 
is also shown. 
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Fig. 6. Directed random market implemented on networks. Numerical simulation considering 
N = 5000 agents and averaging over 250 networks for each type (a) and (b). Every agent starts 
with rrii = (m) = 0.5 and the exchange rule is given by equation LEFT: Stationary probability 
distribution of money P(m), computed after N 2 interactions and entropy-evolution for the whole 
time of the simulation with inset showing transitory regime at the beginning. We can appreciate 
how the system evolves towards a highly entropic state and then it start a self-organising process 
which is analogous to what we found for the gas-like approach in figure [3] but, now, it reaches a 
stable regime. The steady state is characteristic for each network. This comes from the different 
values at which S/S m ax stabilizes. RIGHT: Mean money per economic agent computed according 
to l |17| l and its standard deviation, given by \18\ . Linear fit is obvious. Typical degree distribution 
is also shown. 
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4. Conclusion 

We have found that topology does not determine the final equilibrium distribution 
for the family of undirected random markets, both with or without uniform savings, 
as can be seen by comparing the distributions P(m) from figures [4] and [5] to those 
from figures [I] and [2j Thus, from the uniform value of (m)(k) and er m (fc), we can 
conclude that the connectivity of an agent immersed in undirected random markets 
does not determine if the agent will have more or less money, that is, the degree k 
of an agent does not decide its richness. 



When we consider the directed random market given by the operator ( 14 ) , sim- 
ulations plotted in the figure [6] suggest that this model reproduces very clearly the 
real insight based on the impression that in certain societies the richness is owned by 
a very small fraction of the population and poverty is extended among the majority 
of the agents. For this type of economies, we also discover that the underlying topol- 
ogy determines the stationary distribution of money among the agents, P(m), and 
that the connectivity of each agent proportionally determines its average richness. 

Assuming that we can apply this rule for any quantity that can be exchanged, 
and not only money, this result introduces the interesting idea of how we can create 
systems with a property shared by its inner components with an steady distribu- 
tion essentially determined by the topology defined by the connections between the 
agents, although evidently this statistical equilibrium is dynamical and presents a 
continuous flow between those agents. It is also interesting to highlight how this 
system evolves from the initial zero-entropy state to a state with maximum entropy 
and then how it relaxes towards the asymptotic equilibrium state. 

Let us conclude by saying that we have introduced a new directed random market 
model in the context of economic gas-like models, which can be understood as a 
first and simple model characterized by more connectivity implies more money. 
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